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We derive a continuum theory for the phase transition in a classical dimer model on the cubic 
lattice, observed in recent Monte Carlo simulations. Our derivation relies on the mapping from 
a three-dimensional classical problem to a two-dimensional quantum problem, by which the dimer 
model is related to a model of hard-core bosons on the kagome lattice. The dimer-ordering transition 
becomes a superfluid-Mott insulator quantum phase transition at fractional filling, described by an 
SU(2)-invariant continuum theory. 

PACS numbers: 64.60.Bd, 64.70.Tg, 75.10.Hk 



The standard model of symmetry-breaking phase tran- 
sitions, both classical and quantum, is the Landau- 
Ginzburg- Wilson (LGW) theory Q, where the critical 
properties are described by a continuum theory written 
in terms of the order parameter of the transition. It 
has recently been argued, however, that in certain two- 
dimensional quantum systems, continuous phase tran- 
sitions are possible between symmetry-breaking states 
with apparently unrelated order parameters, in conflict 
with the LGW paradigm 

Another class of non-LGW transitions occurs in clas- 
sical systems with constraints that prevent a fully disor- 
dered state [H- As the temperature is raised, these sys- 
tems instead enter a 'Coulomb phase', where correlation 
functions have power-law forms and strong directional 
dependence. A naive application of the LGW theory 
fails to capture these correlations and so cannot describe 
a transition between an ordered phase and a Coulomb 
phase [||. 

Recent numerical work [HQ indicates that such a tran- 
sition exists in a classical dimer model on the cubic lat- 
tice. This model describes the statistics of close-packed 
'dimers', objects that occupy two neighbouring sites of 
the lattice, with every site of the cubic lattice covered by 
precisely one dimer. There are many configurations that 
obey this constraint, and if all are given equal weight, the 
system displays a Coulomb phase [9( . If instead they are 
given Boltzmann weights that favour parallel dimers, the 
system orders at low temperatures; the ordered phase is 
a six-fold degenerate crystal, breaking the lattice sym- 
metry. 

In this Letter, we outline two steps that lead to a con- 
tinuum description of this classical dimer transition. Our 
first step uses the standard mapping between classical 
statistical mechanics in 3D and quantum mechanics in 
2D, and so provides a bridge between the two classes 
of proposed non-LGW transitions. In the second step, 
we show that long-wavelength properties at the result- 
ing 2D quantum transition are described by the SU(2)- 
symmetric noncompact CP 1 (NCCP 1 ) model. This con- 
clusion is consistent with earlier suggestions, on the ba- 



sis of results for several three-dimensional (3D) quan- 
tum models at finite temperature [H, 0] , that the classical 
dimer transition should be described by a gauge theory 
coupled to multiple matter fields. 

Our approach is to identify the [111] direction as imagi- 
nary time and map to a model of hard-core bosons on the 
kagome lattice. A related mapping has previously been 
applied to another classical transition that is believed 
to lie outside the LGW paradigm, the Kasteleyn transi- 
tion of spin ice in a [100] magnetic field Q. As in that 
case, the Coulomb phase of the classical problem maps 
onto a superfluid phase for the quantum bosons, and the 
long-range power-law form for the correlation functions is 
reproduced by the fluctuations of the Goldstone (phase) 
mode. The dimer crystal becomes a bosonic Mott insu- 
lator, with a chargc-density-wave order that breaks the 
symmetry of the kagome lattice. There is a six-fold de- 
generacy in the ordering pattern, corresponding to the 
degeneracy of the crystalline phase. 

In contrast to the case of spin ice @, the quantum 
phase transition that results is still not amenable to a 
straightforward application of the LGW theory, and is 
instead an example of the class of non-LGW quantum 
transitions described above. The transition is between 
a bosonic superfluid and a Mott insulator at fractional 
filling, and is of the type considered by Balents et al. 0] . 
They showed that, using a dual representation in terms 
of vortices, the continuum behaviour can be described by 
a gauge theory whose form is strongly constrained by the 
symmetries of the original boson problem. In this dual 
model, the condensed phase of bosons becomes the non- 
condensed phase of the vortices, where the fluctuations 
of the photon mode of the gauge theory reproduce the 
superfluid correlations. The condensation of the vortex 
fields leads, by the Anderson-Higgs mechanism, to a gap 
for the photon and short-ranged correlation functions. 

In our case, the continuum theory is an SU(2)-invariant 
gauge theory of a two-component complex vector ip min- 
imally coupled to a noncompact U(l) gauge field A, 

Co = |(V-iA) V | 2 + s |^| 2 + u (|^| 2 ) 2 + K |Vx A| 2 . (1) 
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This can be rewritten in the usual form for NCCP 1 by 
replacing the terms s and u by a hard constraint \tp\ = 1. 
Although the method picks out one direction as imag- 
inary time, space-time isotropy is restored in the final 
form of the theory, a 3D classical continuum model (V 
denotes the three-dimensional derivative operator). 

High-precision Monte Carlo data Q suggest that the 
dimer model has either a single continuous transition as 
the temperature is varied, or a very weak first-order tran- 
sition, with a large but finite correlation length. In both 
cases, we expect there to be a well-defined continuum 
regime with correlation length much longer than the lat- 
tice spacing, described by Eq. ([1]). It remains contro- 
versial whether this continuum theory itself exhibits a 
continuous transition 

An outstanding puzzle noted in Refs. [i| and Q, which 
we are unable to resolve, is that the critical exponent 
values extracted from simulations of the dimer model Q 
differ from those obtained in direct simulations of the 
NCCP 1 model [H [H, [3] and related quantum spin 
models [13, El ■ It has been suggested 0, 0] that the 
discrepancy may be due to the proximity of a tricritical 
point in the phase diagram of a generalized dimer model; 
further numerical work is required to clarify this point as 
well as the nature of the transition. 

The configurations of the classical model are close- 
packed arrangements of hard-core dimers, described 
by the variables d M (r) G {0,1}, giving the occu- 
pation number of the link between the cubic-lattice 
sites r and r + 5^, where S ^ is a unit vector (/i G 
{x, y, z}). The close-packing constraint can be expressed 
as ^ [d M (r) + d M (r — <$ M )] = 1, for all r. The energy of a 
configuration is £ = — nil , where n» is the number of pla- 
quettes (of any orientation) with parallel dimers, though 
our continuum theory applies equally to other potentials 
that also favour columnar crystalline order. 

At zero temperature, T = 0, the dimers maximize n\\ 
by selecting one of the six degenerate columnar ordering 
patterns, distinguished by the staggered order parameter, 

m fl (v) = ±(-iy»[d^r)-d fl (r-5 fl )]. (2) 

The six ground states have m(r) G {±S X , ^S y ,zL6 z }, for 
all r. For T = oo, (m) = 0, and the dimer-dimer corre- 
lation function has the 3D dipole form [§] 
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where ?y r = (— l)^-^ is ±1 on the two sublattices. 

Monte Carlo simulations show a transition, apparently 
continuous, from an ordered phase to a Coulomb phase as 
T is increased, at a critical temperature of T c ~ 1.675 
For T < T c , the lattice symmetry is broken and (m) takes 
on a nonzero value (oriented along one of the cubic unit 
vectors ±<5 M ). For T > T c , the order parameter vanishes 
and the correlation functions have a dipolar form. 




FIG. 1: The projection of the cubic lattice onto a (111) plane, 
with the kagome lattice superimposed. The cubic sites are 
shown by the numbers 0, 1 and 2, giving r mod 3 (see main 
text). Points show the sites of the kagome lattice, situated at 
the centres of the cubic bonds between sites with r mod 3 = 1 
and 2; they therefore lie in planes with r mod 3 = |. The 
larger red circles shown superimposed on some of the kagome 
sites illustrate the occupied sites in one of the six degenerate 
ordering patterns, corresponding to the six ordered states of 
the cubic dimers. The elementary unit vectors of the kagome 
lattice, ei and e2, are shown with dashed blue arrows. The 
coordinate axes shown in the bottom-right of the figure are 
the projection of the cubic x-, y- and 2-axes onto the (111) 
plane. 



We now outline a derivation of the continuum the- 
ory for this transition, based on a mapping to quantum 
bosons in 2D. Our first step is to define the imaginary 
time t = r^, given by the projection of the 3D posi- 
tion onto the [111] direction. We then map the configu- 
rations of a given (111) plane of the classical problem to 
the basis states in a quantum Hilbert space, by simply 
identifying the presence (or absence) of a dimer with the 
presence (absence) of a boson. As illustrated in Figure [IJ 
the midpoints of the cubic bonds satisfying t mod 3 = § 
form a set of kagome planes stacked in the [111] direc- 
tion, and we use this set as the space-time lattice for the 
quantum problem. 

For this mapping to make sense, we require conser- 
vation of particle number. The close-packing constraint 
implies that 



EE^)= E 



re(r-l) L 



l-E^(r) 



(4) 
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where r £ r indicates a sum over all cubic sites in 
imaginary-time slice r. This implies that if the kagomc 
plane at r = | has n bosons and a total of N sites, then 
the plane at r = | + 3 will have ^ — n bosons. We must 
therefore choose an imaginary-time step St that is even 
in order to give particle conservation; the smallest choice 
is St = 6 and so we take only the bonds whose midpoints 
are at r mod 6 = § to define the quantum problem. 

The six distinct ground states of the dimers map to 
six ordered states of the kagome bosons, and the exam- 
ple with m = — S z is shown in Figure [1] In each ordered 
state, half of the sites on a given sublattice are occu- 
pied, so the density (bosons per site) is / = g. The low- 
temperature phase of the classical model maps to the or- 
dered phase of the bosons, while in the high-temperature 
phase the bosons condense and there is no spatial order, 
corresponding to the Coulomb phase with (m) = 0. 

The basis states for the quantum Hilbert space are 
given by the occupation-number states of hard-core 
bosons on the kagome lattice. The hard-core nature of 
the dimers further implies that there is hard-core nearest- 
neighbour repulsion, and each triangle of the kagome lat- 
tice (of either orientation) can be occupied by at most one 
boson. The partition function for the classical problem is 
Z = TtT l / St , where L is the length of the system in the 
[111] direction and periodic boundary conditions are as- 
sumed in this direction. The transfer matrix T has rows 
and columns labelled by the configurations of two planes 
separated by St = 6, and its elements give the statistical 
weights for these configurations, summed over all possible 
configurations of the intermediate planes. The effective 
quantum Hamiltonian Tt is defined by T = e~ nSr , so 
that Z is given by the quantum partition function at in- 
verse temperature (3 cx L. The classical thermodynamic 
limit is therefore given by the quantum zero-temperature 
limit, f3 — > oo. 

It is in principle possible to find T, and hence H, ex- 
actly for a finite lattice, by enumerating all possible con- 
figurations of the dimer problem. For even quite small 
lattices, however, this is a computationally difficult prob- 
lem, and we have made no attempt to find T explicitly. 
Instead, the continuum theory for the transition can be 
found by using general principles such as symmetry to 
place constraints on 7i. 

The Hamiltonian will take the form of a generic lattice- 
boson model, with all hopping terms and interactions 
that are consistent with symmetry, as in the related map- 
ping for spin ice Q . The parameters in the Hamiltonian 
are functions of those in the classical model, and so in our 
case the hopping and interaction strengths will depend on 
the classical temperature T. Since the high-temperature 
phase corresponds to a boson condensate, we expect the 
ratio of the hopping strength to the interactions to in- 
crease with T. 

For the mapping to bosons to be useful in the thermo- 
dynamic limit, we require that the Hamiltonian should be 



local, at least in its action on low-energy states. We have 
no general proof that this is the case, but we argue that it 
is so in the region of interest, near the transition. In this 
case, the low energy states of the bosons can be written 
in terms of states with a finite density of domain walls 
separating regions with different density-wave orderings. 
These one-dimensional domain walls are, in 3D, the in- 
tersections of a given (111) plane with two-dimensional 
domain walls between dimer orderings. Consider a do- 
main wall that moves by a large amount in a single time 
step, so that one of the two neighbouring domains grows 
by an area SA. In the three-dimensional picture, this do- 
main wall has a section of area SA that runs parallel to 
the (111) plane. Such a configuration costs entropy and 
configuration energy that grows with SA, and hence has 
an exponentially suppressed contribution to the transfer 
matrix. 

We defer a full analysis of the symmetry properties of 
the Hamiltonian TL to a forthcoming paper, and discuss 
here only a symmetry of particular importance: rotation 
by | about a [111] axis passing through the centre of 
a kagomc hexagon, and hence through a cubic lattice 
site with r mod 3 = (such as the one at the centre of 
Figure [J). This maps the kagomc lattice onto itself, but 
is not a symmetry of the cubic lattice, since it exchanges 
the positions of the cubic sites on the intermediate planes. 
(For instance, it exchanges cubic sites with r = 1 and 
r = 2.) We therefore define the operation R consisting of 
this rotation followed by a reflection in the (111) plane 
r = |, which is a symmetry of the cubic lattice. In terms 
of the bosons, M is a rotation followed by a time-reversal 
operation t — > 3 — r, and it commutes with 7i. 

Before addressing the transition, we note that tak- 
ing full account of the symmetries allows an alternative 
derivation Q, starting from the quantum model, of the 
dipolar form of the Coulomb-phase correlation functions. 
This phase corresponds to the superfluid, where long- 
range correlation functions are dominated by fluctuations 
of the Goldstone phase mode 0, described by the coarse- 
grained action £^ ~ |c?0| 2 . Correlation functions can 
be found by using the symmetries to express the boson 
number operators and dimer occupations in terms of the 
continuum field 0, giving for instance ~ rj r d^<j). 

The continuum theory describing the transition be- 
tween a superfluid and a Mott insulator at fractional fill- 
ings can be expressed in terms of dual vortex fields. We 
will sketch the derivation of this theory; readers are re- 
ferred to the papers of Balents et al. 0] for an overview of 
the approach. In essence, it consists of a transformation 
from a current-loop representation of the boson problem 
[lot , where the basic degrees of freedom are the currents 
Jn(r) defined on the links of the space-time lattice, to a 
gauge field A^(r) defined on the links of the dual lattice, 
according to J = V x A. 

Note that the bosons belong on the sites of the kagomc 
lattice, and so the space-time lattice is not the original 
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cubic lattice of the dimer problem, but instead consists 
of stacked kagome planes. The dual of kagome is the dice 
lattice [HI [lij , and so the gauge field A is defined on the 
links of a lattice of stacked dice planes. 

The dual theory has a gauge invariance resulting from 
the redundancy of the definition of A, and can be writ- 
ten in terms of the gauge field ^4 M (r) and matter fields 
?/>(r) defined on the sites of the dual lattice. In terms 
of the original bosons, these matter fields represent vor- 
tices, and the physical density of bosons / = ^ becomes 
an effective magnetic field Vxi for the vortices. The 
low-energy properties of the theory are then determined 
by the vortex dynamics in this effective magnetic field. 

For bosons at filling / = g on the kagome lattice, the 
conclusion of this analysis is that there are two degener- 
ate minima of the vortex dispersion within the reduced 
Brillouin zone The continuum theory can then be 

written in terms of field operators tpo and tpi correspond- 
ing to these low-energy modes, and the symmetries again 
strongly constrain the possible terms in the action. 

An important example is the symmetry under R, which 
transforms the fields according to 

(<f>0\ RC^fl 1\ fip \ 

implying that this must be a symmetry of the continuum 
theory. While the full action is only symmetric under a 
discrete subgroup of SU(2), the symmetry under R for- 
bids any terms that break SU(2) up to (and including) 
fourth order in the fields <p. 

The allowed interaction terms can be found by express- 
ing all gauge-invariant bilincars of the fields ip and ipi in 
terms of the boson density operators Q , or equivalently, 
in terms of the dimer order parameters m M . Using the 
symmetries of the cubic lattice, we find 

~ (p^cr^cp , (6) 

where cr^ are the Pauli matrices. 

Defining the SU(2) Casimir invariant £1 = l^oP + l^il 2 ; 
all gauge-invariant combinations of the fields <p can be 
written in terms of Q and m x y z , and it is straightforward 
to show that fl 2 ~ |m| 2 . While the action can contain 
any term involving only £1. terms involving functions of 
m M that break SU(2) are strongly restricted by symme- 
try. The first such term allowed by Eq. |[5j), of sixth order 
in (p, is m x m v m z , and is excluded by requiring symmetry 
under m — > — m. The lowest-order combination satisfy- 
ing all the symmetries is of eighth order, and is given by 
C\ = v m^j. In order to describe the transition to the 
particular ordered states with which we are concerned, 
we require v < 0. 

The SU(2)-invariant part of the action, given in 
Eq. (fT]), takes the standard form for a complex vector 
minimally coupled to a noncompact U(l) gauge field. 
We have omitted terms containing higher derivatives or 



higher powers of the field ip. The cubic symmetry of the 
original dimer problem means Co must be space-time 
symmetric, and we have made this explicit in Eq. ([1]). 

It is not firmly established whether this action has any 
nontrivial fixed points under the RG. As noted in Rcf. 3, 
a conjecture of this sort is difficult to test analytically 
[ill, and the evidence from numerical studies is still in- 
conclusive. It is clear from simulations, however, that 
the transition is, at most, weakly first order, with a very 
large correlation length, and so a continuum description 
is appropriate. The fact that L\ is of eighth order in the 
field ip makes it highly likely that this term is irrelevant 
in the continuum, and that the effective theory has an 
emergent SU(2) = 0(3) symmetry. 

Using this assumption, and the expression for the crys- 
talline order parameter given in Eq. (|6|), we are able 
to explain two qualitative observations made by Mis- 
guich et al. Q based on their numerical results near the 
transition. Firstly, the data strongly suggest that the 
order-parameter distribution is spherical, as expected for 
an 0(3)-symmetric theory. Secondly, the dimer-dimer 
correlations are dominated by a 'spin-like' contribution 
This follows from the fact that (and hence the 
dimer variables g? m ) couples directly (without derivatives) 
to a bilinear in the critical field ip, and forms a three- 
dimensional representation of SU(2). We therefore find 
(m^m v ) ~ S fl4/ \r\^ d+2 ^ Tlm The absence of a weak- 
coupling fixed point again prevents us from making quan- 
titative predictions about the anomalous dimension rj m . 
As noted by Misguich ct al. [||, these properties, while 
explained straightforwardly by the NCCP 1 theory, are 
incompatible with other, more obvious, candidate con- 
tinuum theories, such as the 0(3) model. 

In conclusion, we have derived an NCCP 1 theory to de- 
scribe the observed phase transition in a classical dimer 
model on the cubic lattice 0], by means of a mapping 
to a quantum transition in two dimensions. Simple qual- 
itative predictions based on this mapping are in agree- 
ment with numerical results for the dimer system [5], [|| . 
The discrepancies between the observed properties of the 
NCCP 1 and dimer models, as well as the uncertainty re- 
garding the nature of the transitions, illustrate the need 
for further numerical studies of both models. 
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